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Abstract. We consider a quantum state interacting with a short intense linearly 
polarized laser pulse. Using the two-dimensional time representation and Floquet 
picture we establish a straightforward connection between the laser carrier-envelope 
phase (CEP) and the wave function. This connection is revealed as a unitary 
transformation in the space of Floquet components. It allows any CEP effect to be 
interpreted as an interference between the components and to put limits on using the 
CEP in coherent control. A 2-level system is used to illustrate the theory. On this 
example we demonstrate strong intensity sensitivity of the CEP effects and predict an 
effect for pulses much longer than the oscillation period of the carrier. 



1. Introduction 

Progress in manipulating ultrashort intense laser pulses has made possible studying 
laser-matter interactions in qualitatively new regimes. For instance, very short pulses 
having only a few oscillations of the laser field can be produced |J] 12] . In contrast to the 
more conventional case, when a laser pulse is much longer than the carrier period, 
the carrier phase of the pulse with respect to the envelope maximum can become 
an important parameter for short pulses. This phase is called the carrier-envelope 
(CE) phase (CEP). It has been demonstrated experimentally, that the CE phase can 
significantly influence ionization of Kr atoms by infrared laser pulses [TJ. Recently, 
similar experiments were performed with Rydberg states of Rb atoms ionized by a few- 
cycle 25 MHz pulse [3], and the spatial distribution of the ionized electrons has shown 
sensitivity to the CEP. Even potentially stronger effects were predicted theoretically 
for dissociation of the HD + in the laser field [4] and experiments are being performed. 
The sensitivity of high harmonic generation (HHG) on the laser CEP is also known 
[2|. Molecular isomerisation in short intense laser pulses also provides an interesting 
example where CEP effects are important [5|. 

However, the CEP effects are, probably, not fully understood theoretically. Only a 
few models that go beyond the qualitative picture have been discussed. An interesting 
interpretation of the CEP effects in ionization as a double-slit interference in the time 
domain was proposed by F.Linder et al. [6]. This interpretation, however, does not 
help to describe results for high energy electrons, or, more generally, to describe the 
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dependence of the CEP effects on the final state energy. The goal of this paper is to 
propose a simple and physically significant interpretation for the influence of the CE 
phase on the interacting wave function. We propose a picture revealing a simple but 
exact relationship between the laser phase and the state. 

We use the 2D time (or (t,t')) formalism [7] together with a Floquet representation 
of the interaction. This approach allows a model-free separation of different time scales 
present in a short laser pulse. We shall see how the CEP can be eliminated from the 
equations by a simple unitary transformation. This transformation allows recovery of 
the CEP dependence of a final state from the wave function components for one CE 
phase only. 

We illustrate our result by demonstrating how the state of a two-level system (a 
qubit) can be controlled by a short pulse, revealing the mechanisms of CEP influence 
on the state of the system. We demonstrate and explain dependence of the CEP effect 
magnitude on the maximal field of the pulse and the pulse duration. It is shown that 
CEP effects can be observed even for pulses that are much longer than one oscillation 
period of the carrier. 

2. Theory 

We base our approach on investigating the time-dependent Schrodinger equation 

i^!t = [H Q + V(t)}9. (1) 

Here \I> is the wave function, H is the Hamiltonian of the system without laser field, and 
V(t) stands for the laser-matter interaction potential. We shall consider laser-matter 
interactions of the following form 

V(t) = -E(t) ■ d cos{ut + (p) . (2) 

Here E(t) is the envelope of the laser pulse field, d is the dipole interaction operator, 
uj is the laser carrier frequency and p is the carrier-envelope phase. The latter can be 
a very important parameter especially for sufficiently short pulses, and it is the main 
parameter studied in this work. 

In what follows, we give a brief description of the 2D time formalism [7] and 
introduce a two-dimensional time representation for a system in a pulsed laser field. 
We shall demonstrate how this representation allows eliminating of the CEP p 
from the evolution equation and construction of a CEP-independent solution. The 
CEP dependent solution of the initial equation @ can be recovered from the CEP- 
independent solution with a unitary transform. 

2.1. The 2D time formalism for a system in a periodic external field 

The formalism of two-dimensional time is very useful to treat time-dependent systems 
that show both periodic and non-periodic behavior, such as atoms and molecules in a 
field of a laser pulse. One of the advantages of this formalism is that it allows separation 
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of periodic and non-periodic dynamics without resorting to adiabatic expansions that 
might converge slowly for nonadiabatic systems. 

In the 2D-time representation one introduces a second time s, such that the envelope 
and the periodic factors in the laser-matter interaction (2) depend on different time 
coordinates 

V(s,t) = -(E(s),d)coa{ut + <p) 

H(s,t) = H + V(s,t) . [ ' 

One also introduces a second time dynamics 

(i— + i-Ms,t) = (Ho + V(s,t;^Ms,t) . (4) 

It is not difficult to see that if \I/(s, t) satisfies (j3j), the solution of (DO) can be written 
as *ff(t) = *&(s,t) \ s=t . Indeed, restricting the solution to the "diagonal" time s = t and 
substituting it to the left-hand side of the equation (TJ) we get the left-hand side of the 
equation (3} 

The right-hand sides of the equations (GO) and (3J are identical at s = t. Thus, once 
the 2D-time equation is solved, one has a solution of the original equation (1) as well. 
More detailed discussion of the 2D-time formalism can be found at the original paper 
of Peskin and Moiseev [7J. 

2.2. Floquet representation for a finite-time pulse 

If the laser pulse duration is not considerably shorter than the oscillation period 2n/u, 
expanding the wave function \l/(s, £) into a Fourier series in t can be a reasonable way 
of solving equation (j4j, even for pulse duration comparable with the oscillation period 
\E\. In fact, this approach holds as far as the description of the laser pulse itself in terms 
of carrier and envelope is valid pTj. In this paper we discuss pulses longer than one 
oscillation period, and, thus, our approach to solving equation (0| is justified. 

We start constructing the Floquet representation [9J, TO] for the laser pulse from 
expanding the wave function 

oo 

*(M)= ^n{s) ■ (5) 

n=~ oo 

We shall call the coefficients 4> n (s) n-photon emission (n > 0) and absorption (n < 0) 
amplitudes. 

Bringing the time derivative zj^ to the right hand side of equation (3J and 
substituting the representation @ we come up with the following infinite system of 
equations 

ifM s )= ie*(E(s),d)^_i(s) 
+(H + nu) (j) n {s) 
+|e-^(E( S ),d) ct> n+1 {s) U 
n = — oo, . . . , — 1, 0, 1, . . . , oo 
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Starting at time s = s_ from the initial state 

(po(s-) = ip 
Ms-) = n ± 

and propagating the amplitudes according to (6) sufficiently long, we can investigate 
different laser-induced processes, such as dissociation, ionization, (de) excitation etc. 

It is convenient to introduce a vector of n-photon amplitudes <&(s; </?) and a Floquet 
Hamiltonian 

H( S ;^) = Ho + V(^) 

where H stands for the diagonal and V(s; <p) for the off-diagonal terms in the equation 
(MI). The equation @ can be written as 
d 

i— ®{s;ip) = H(s;y?)$(s;y?) . 

In this notation we emphasize that the Floquet Hamiltonian and the solution both 
depend on the carrier-envelope phase (p. This dependence, however, is parametric, and 
it can be eliminated. 



2.3. Unitary equivalence of the CEP 

The obvious advantage of equation {6j for studying the CEP effect is that the phase 
dependence enters the equation linearly. As we shall see, it allows exclusion of the CE 
phase from the Floquet Hamiltonian by a simple unitary transformation. 

Let us introduce the following operator U((p) acting in the n-photon amplitude 

space 

[l%)] mn = 8 mn e- in * . 

Now consider the Floquet Hamiltonian H(0)corresponding to <p = 0. It is easy to verify 
that the following relation holds: 

Hfa) = Utfe>)H(0)UfoO 

and to rewrite the evolution equation as 

^ S U(^)0( S ;^) = H(O)U(^)0( S ;^) . 

This equation helps to establish the main result of this paper, which is the unitary 
equivalence of solutions corresponding to different carrier-envelope phases 

*(s ;¥ = U%)*(3;O) ■ (7) 

This relationship allows any CEP effect to be interpreted as interference of different n- 
photon channels. It can be clearly seen if we recall the expression for the wave function 
in physical time 

oo 

E e *"^n(*;0) , (8) 
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where we have introduced the final state components of the wave function 4> n (t; 0) = 
e inuJt (p n {t;0). After the laser pulse is off at the moment t = t , the components 
4> n depend on time uniformly due to the internal dynamics of the system only, i.e. 
<f)(t) = e- iHo(t - to) 0Oo)- Depending on the laser CEP, the wave function components (p n 
gain phases proportional to the corresponding net number of exchanged photons n and 
the CEP p. 

It is important to note that all the given results are obtained under rather general 
assumptions of dipole laser-matter interaction and a stable shape of the laser pulse. We 
did not assume anything specific for any particular physical system. That means that 
any CEP effect can be considered as interference of several n-photon channels. 

2.4- CEP effect observation 

Before discussing application of the formula f8j to particular model systems, we say a 
few words on observing the phase effects in general. 

Let O be an observable of interest. Using the representation (8) at large times we get 
an explicit CEP dependence of the mean value of O through the n-photon components 
at zero phase 

oo 

(6)= e^-^(0 n (t o ;O)|O(t)|0 fc (t o ;O)) , 

k,n=— oo 

where 0{t) = e in o(t-to)Q e -iUo(t-to) ^ Rearranging the terms in the series and using 
& = O we get a Fourier series for the (O) CEP dependence 

oo 

(O) (p) = — + ^^(Re«fc cos kp + Ima^ sin kp) (9) 

k=l 

with the coefficients defined as 

oo 

a k = 2 (<f>n-k(t ] 0)\6(t)\(f> n (t ; 0)) . 

n=— oo 

It is useful to introduce a measure of CEP effect observability. This measure can 
be chosen as a norm of the (^-dependent part in equation Q 

oo 

o = (£>*f)* • (10) 
fc=l 

We shall refer to this quantity as absolute CEP amplitude of the observable O, since it 
indicates how much (O) can deviate from its mean value «o when varying the CEP. It 
is worth mentioning that a is directly connected to the mean square deviation of (O) 
from its CEP-averaged value, namely: 
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It is also useful to have a CEP amplitude weighted with the mean value a and 
experimental sensitivity 5 

3. Numerical illustrations 

In the following section we give two types of demonstrations: quantitative and 
qualitative. 

Quantitatively, we shall check the agreement of equations (Q| and @ with CEP 
results calculated independently. For that type of demonstration we have to calculate 
the n-photon amplitudes either from a direct solution of equation (4) or by extracting 
them from an independently calculated final state wave function. Knowledge of the 
amplitudes, however, is rather expensive. To predict the CEP response of any particular 
physical system, we have to know not only the population probabilities of the Floquet 
states, but their phases as well. Usually, this information can be obtained only from a 
numerical solution of the TDSE, and we shall explain how expression (j8| can be used 
to reduce the amount of numerical calculations. 

Qualitatively, we shall discuss the physical conditions needed to observe CEP effects 
experimentally. Such conditions are equivalent to the existence of several interfering 
components. In fact, as equations (HJS) suggest, CEP effects exist if and only if several 
n-photon components contribute to the same physical state. We shall see how this 
condition is realized in a simple 2-level model. On this example we shall discuss intensity 
and pulse duration dependence of the CEP effects. We should mention, however, 
that not all the physical systems conventionally treated in a 2-level approximation are 
suitable for clear CEP effect demonstrations. For instance, choosing an experimental 
realization, one has to fulfill the condition that the two states used in the demonstration 
must be well separated from other eigenstates of the system, because the important 
transitions have essentially nonresonant multiphoton character. Such systems can be 
realized, for example, as double quantum dots [US, [T4] or as ionic hyperfine qubits 



15] . Although the 2-level model cannot describe processes of ionization or molecular 
dissociation realistically, it still gives a reasonable qualitative description of the CEP 
effect observability. More detailed study of the CEP effects involving continuum states 
is a subject for another investigation. 

In all the examples we shall use a Gaussian shape of the laser pulse 

V(t) = -d-E(t) = V Q e~ { r )2 , 



where E is the peak field, \/2 ln2r is the intensity FWHM pulse duration and 
Vo = —d ■ E(T) is the peak interaction energy. Since we are particularly interested 
in molecular systems, we choose energy scales in the typical energy range of molecular 
vibrational states and fix the laser carrier frequency to 0.058 a.u., what corresponds to 
a standard 790 nm Ti:Sapphire laser. 
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3.1. Excitations in a 2- state model. 



We start from a simplest example of an excitation in a 2-state system. Let \a) and \b) 
be the two eigenstates with energies E a and E\,. Without loss of generality we can set 
the first energy to zero, E a = 0, Eh = E a + A = A, and consider \a) the initial state. 
The wave function in this case takes form 

il> = a{t;<p)\a)+b(t;(p)\b) . 

Suppose the coupling to the laser field between the two states is proportional to the 
laser field, such that the corresponding time-dependent Scrodinger equation reads 



.d_ ( a(t) 
l di [ b(t) 





V(t) COs(iUt + (f) 



V(t) cos(uot 
A 




.(12) 



After introducing the coupling energy V(t) = —d - E(t) and keeping only states coupled 
to the initial state \d) the Floquet Hamiltonian H for p = in (H) can be written as 



H 







V(s) 
2 





















V(s) 
2 


- 2uo 


V(s) 
2 





















V(s) 
2 


A-u 


V(s) 
2 





















V(s) 
2 





V(s) 
2 





















V(s) 
2 


A + uj 


V(s) 
2 





















V(s) 
2 


+ 2u 


V(s) 
2 




V 

















V(s) 
2 




J 



In this representation even Floquet amplitudes of correspond to the ground state \a) 
and odd ones to the excited state \b), i.e. 

|02n) = a 2n (t)e- 2n ^\a) , 
|02n+i) =b 2n+1 (t)e-( 2n+1 ^\b) , 
n = — oo, . . . , oo . 

Let us demonstrate first how the CEP reveals itself in the excitation probability. At 
first, we calculate the Floquet amplitudes for ip = numerically. In this example we use 
the energy gap A = 0.066 a.u., which slightly exceeds the photon energy, and the dipole 
coupling is chosen as d = —1 a.u. The pulse intensity is chosen as 5 x 10 14 W/cm 2 , which 
corresponds to the peak interaction energy V = 0.1194 a.u., and the pulse duration is 
t = 248 a.u., which is about 7 fs intensity FWHM. The calculations were performed 
with 30 Floquet blocks, what guarantees accurate results even for intensities higher 
than 10 15 W/cm 2 . Let us look for the excitation probability P ex as the observable. 
We choose the final propagation time t = 4r = 992 a.u., when the field is negligible. 
According to equation (BJ), the final state wave function is expressed in terms of the 
Floquet amplitudes as 



+ 



Eoo 
n- 



3 i2n ^ 2n (t )) | a) 
e« 2n+1 H 2n+1 (t )) e- A ^\b) 



(13) 
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Figure 1. 2-level system photo-excitation probability as a function of the CEP for the 
energy gap A = 0.066 a.u., laser frequency w — 0.058 a.u., peak field Vq = 0.1194 a.u., 
pulse duration r = 248 a.u.. The dots are obtained from the direct solution of equation 
(ECS), the line is constructed from the final state n-photon amplitudes for tp = 0. 

Thus, in order to calculate the excitation probability we have to evaluate the odd n- 
photon components of the final state 

oo 

P ex = \(b\V)\ 2 = I e* 2n+1 H 2n+1 (t )\ 2 

n=— oo 

The following final state amplitudes will contribute to the excited state of the system at 
the final time &i(4r; 0) = 0.0092 - 0.0202i, Li(4r; 0) = 0.1445 - 0.3178i, L 3 (4r; 0) = 
—0.0543 + 0.1193i. All other amplitudes are negligibly small. According to equation 
(1131, the excited state component CEP dependence reads 

b{t- cp)\b) = (b x e^ + L ie - iv + b^ 3 e- 3 ^)e- iA{t - tn) \b) . 

This results in the following explicit expression for the excitation probability: 

P ex = 0.1395 - 0.07602 cos 2y? - 0.00582 cos4y? . 

This line is shown in the Fig. CD together with results of direct solution of equation 
(CCS). Since no approximations were made, besides cutting off the series (CT3l), the perfect 
agreement is not surprising. 

Giving this example, we calculated the Floquet amplitudes by direct solution of the 
system of equations (6). In practical applications, however, this approach is not effective: 
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usually, solving a small system n times is simpler than solving an n times bigger system 
of equations once. In the case of the equations @ the situation is even worse. Even if 
we have only 3 components contributing to the final state, as in our example, one needs 
more than 5 Floquet blocks kept in the equations to reproduce the correct dynamics even 
for moderate peak fields. Because of that, recovering the amplitudes by subsequently 
solving the original Schrodinger equations directly for several CE phases and fitting 
the expression @ to the wave function should be a much more effective procedure for 
realistic calculations. 

After giving this numerical example, we are ready to discuss some qualitative 
properties of the CEP effects. Are there any conditions for the peak interaction energy 
Vo and the pulse duration r that limit an observability of CEP effects? 

Let us discuss the intensity dependence first. In order to have interference, one 
has to make sure that there are several n-photon components contributing to the states 
with the same final state energy. It cannot be achieved with one-photon transitions 
only, so there must be a minimal intensity that allows a clear observation of the CEP 
effect. To understand the lower intensity limit, consider the eigenvalues of H shown 
in Fig. 121 On the leading edge of the laser pulse the interaction energy grows up, this 
correponds to moving from the left to the right in Fig. El After the peak, which defines 
the rightmost point in Fig. \2{ the system goes from the right to the left during the 
trailing edge of the pulse. As one can see from tracking the eigenvalues of the dressed 
system, the two crossings that we need to populate the state \b) via 1- and 3-photon 
transitions correspond to the interaction energies Vti ~ a.u. and Vt2 ~ 0.1 a.u.. 
This gives us a limitation on the laser peak intensity: no interference can happen if the 
maximal field-matter interaction energy does not approach the second crossing. We can 
safely say that in our example no CEP effect is expected if the peak interaction energy 
is smaller than 0.05 a.u.. This observation is illustrated in Fig. \E\, where we have plotted 
the weighted CEP effect amplitude (TTI). As we expected, we see an essential growth 
of the CEP effect contrast only above 0.05 a.u., when the peak field starts approaching 
the 3-photon crossing. 

The big maximal intensity alone, however, is not enough to observe the effect: 
the field should change fast enough when passing through both crossings, otherwise 
transition happens adiabatically, and no population transfer occurs. The proper timing 
conditions must be satisfied. This leads us to the question of the pulse duration 
dependence of the CEP effects. 

As we have seen, nonadiabatic second transition (see Fig. [2} is a necessary condition 
for the CEP observation. We can use the Landau-Zener model to qualify the presence 
of the CEP effects by estimating the second transition probability as a function of the 
pulse duration. Let Vr2 be the position of the second crossing. For the energy gap of 
our example (A = 0.066 a.u.) the level splitting is about g(Vr2) = 0.016 a.u4 The 

| For small energy gaps compared to the photon energy we can use a Bessel approximation to the 
spectrum of the Floquet Hamiltonian. This allows us to estimate the splitting of the Floquet eigenstates 




V, a.u. 

Figure 2. 2-level system Floquet Hamiltonian eigenstates as a function of the laser- 
matter interaction energy V = —E ■ d. The model parameters are uj = 0.058a.u., 
A = 0.066a.u. Thick lines correspond to the ground state and dotted lines to the 
excited state. 

growth rate of the interaction energy when crossing the second transition region 
dt r V V T2 

plays the role of velocity in the Landau-Zener formula. Now we are ready to estimate 
the second transition probability as 

g = g(V T 2) 2 r (14) 

2V T 2y/HV /V T 2) 

If no higher order transitions happen, the amplitude of the CEP effect should be 
proportional to the square root of this probability, since it is proportional to the 
amplitude rather than the probability of the m-photon state population. This square 
root of the probability is shown in Fig. \4k together with the amplitude of the excitation 
probability CEP dependence. It is clear that probability ( 1141 ) should have a substantial 
value to make possible an observation of the CEP effects. Equation ( JT4T ) shows that 
CEP effects should disappear exponentially with growing pulse length. 

at the crossing point as g{V T 2) « - #f v/|)| ~ ~ 0.40-D)]. 
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Peak intensity, W/cm 




Figure 3. The weighted excitation probability CEP amplitude aw ([TTT) as a function 
of the peak interaction energy. No phase dependence below Vq = 0.05 a.u. is observed. 
Calculations were performed with 30 Floquet blocks, the sensitivity parameter 6 in 
(CLTj) is set to 0.01. 



It is important to mention that the pulse duration needed to observe the CEP effects 
is a property of the laser-matter interaction rather than the laser pulse alone. In the 
two-level model, as equation ( JT4T ) suggests, it is defined by level splitting at the second 
crossing point and the position of the second crossing. For large energy gaps it is easy 
to obtain smaller multiphoton splitting. This allows one to predict the CEP effects even 
for pulses that are substantially longer than one oscillation period, as demonstrated 
in Fig. 3b. There we show the CEP observability a together with its Landau-Zener 
estimation. Even for pulses longer than r = 1700 a.u., what is about 30 periods in the 
field FWHM, we still can see a variation of the probability with the CEP about 10%. 
We must mention, however, that experimental observability of the long-pulse effect can 
be limited by the stability of the pulse shape, which might be difficult to keep at a time 
scale smaller than one oscillation period for long pulses. 
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Figure 4. CEP variation of the excitation probability as a function of the pulse 
duration. The black lines correspond to the Landau-Zener estimation of the effect via 
the transition probability at the second crossing, a) A = 0.066 a.u. b) A = 0.12 a.u. 
The peak interaction energy is taken about 30% higher than the second transition 
energy. 

4. Summary 

We have investigated a quantum state experiencing a laser pulse of a stable envelope 
shape and varying CEP of the oscillatory part. We have shown how the CEP can be 
excluded from evolution equations and how the CEP dependence of the final state can 
be recovered from CEP-independent results. On the example of a 2-level system we 
have demonstrated a critical intensity and pulse duration dependence of CEP effects. 
In contrast to the common conception that CEP effects can be expected only when the 
pulse duration is nearly as short as the laser oscillation period, we have demonstrated 
that CEP effects can exist even for pulses much longer than that. The pulse duration 
that allows the CEP effect observation critically depends on the properties of the 
system interacting with a laser pulse. In the 2-level system, long-pulse CEP effects 
can be observed in essentially above-threshold excitation regimes, when the energy gap 
is considerably larger than the photon energy. 

The approach which is suggested in this paper is rather general, and it would 
be interesting to study more complex systems from this point of view. For instance, 
it would be interesting to study how n-photon component interference affects high- 
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harmonic generation, ionization and dissociation in different systems. There are also 
many physical and mathematical questions to be addressed, such as gauge-invariant 
formulation of the theory and checking the classical limits of the interference effect. 
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